Level statistics in arithmetical and pseudo-arithmetical chaos 
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We resolve a long-standing riddle in quantum chaos, posed by certain fully chaotic billiards with 
constant negative curvature whose periodic orbits are highly degenerate in length. Depending on the 
boundary conditions for the quantum wave functions, the energy spectra either have uncorrelated 
levels usually associated with classical integrability or conform to the "universal" Wigner-Dyson type 
although the classical dynamics in both cases is the same. The resolution turns out surprisingly 
simple. The Maslov indices of orbits within multiplets of degenerate length either yield equal 
phases for the respective Feynman amplitudes (and thus Poissonian level statistics) or give rise 
to amplitudes with uncorrelated phases (leading to Wigner-Dyson level correlations). The recent 
semiclassical explanation of spectral universality in quantum chaos is thus extended to the latter 
case of "pseudo-arithmetical" chaos. 
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Introduction: After more than two decades of inves- 
tigations, the famous BGS conjecture [l[ has recently 
found a semiclassical explanation The quantum 

level spectra of classically chaotic systems display fluctu- 
ations conforming to the random-matrix-theory (RMT) 
predictions for the Wigner-Dyson universality classes. 

However, there is the notable exception of "arithmeti- 
cal" systems which are fully chaotic classically but dis- 
play quantum spectral statistics close to Poissonian, a be- 
havior usually associated with integrable classical motion 
• Quantum mechanically these exceptional dynam- 
ics exhibit an infinite number of the so called Hecke op- 
erators commuting with the Hamiltonian. Therefore the 
energy spectrum falls into practically independent mul- 
tiplets such that nearby levels bear no correlation. On 
the classical side the periodic-orbit action spectra of such 
systems are distinguished by a degeneracy exponentially 
growing with the orbit period. 

On the other hand, by merely changing the boundary 
conditions for the quantum wave functions of some such 
exceptional arithmetical systems one can retrieve univer- 
sal spectral fluctuations a la Wigner and Dyson while not 
at all changing the classical dynamics. It is customary to 
speak of pseudo-arithmetical systems then 1^- 15j|. 

The strikingly different quantum behavior of arith- 
metical and pseudo-arithmetical systems might raise 
doubts about the validity of the recent semiclassical ex- 
planation of universal quantum spectral fluctuations un- 
der conditions of classical chaos; after all, the classical 
dynamics are identical for the systems under discussion, 
and so appear, on first sight, the Gutzwiller type semi- 
classical periodic-orbit expansions. Various suggestions 
were ventured for the effect of the boundary conditions, 
among them a distinction between the orbit classes con- 
tributing to the Selberg trace formula applicable (and 
exact) in the arithmetical case, and the Gutzwiller for- 
mula applicable in the pseudo-arithmetical case. 

We show here that the explanation is much simpler and 



lies in special properties of periodic orbits. Due to these 
peculiarities all equal-length orbits (save for a negligible 
fraction) of an arithmetical system contribute Feynman 
amplitudes with the same Maslov phase; their construc- 
tive interference makes for nonuniversal spectral statis- 
tics. In the pseudo-arithmetical case these phases vary 
randomly within a degenerate-action multiplct such that 
destructive interference makes the high action degener- 
acy ineffective. The difference between the two cases is 
most easily revealed for the diagonal approximation to 
the spectral form factor, and therefore that approxima- 
tion will play a central role here; off-diagonal corrections 
will be discussed briefly in the end. 

The billiard T*(2, 3, 8): We shall not deal with the al- 
ternative arithmetical/pseudo-arithmctical in full gener- 
ality but prefer to work with a representative example, 
the so called triangular billiard T*(2,3,8). That system 
was first considered in the studies of the free motion on 
the surface of constant negative curvature tesselated by 
regular octagons [lij. Desymmetrization of the regular 
octagon necessary to get rid of the rotational and re- 
flection symmetry leads to a triangular fundamental do- 
main with the angles tt/2, 7t/3, tt/8. Depicted inside the 
Poincare disk \z\ < 1 in the complex plane, the triangle 
has its hypotenuse (N) and the longer leg (L) directed 
along two diameters whereas the shorter leg (M) looks 
like an arc (Fig. 1). A classical periodic orbit folded 
into the triangle looks like a sequence of arcs mirror- 
reflected from the sides of the triangle. The classical mo- 
tion is completely chaotic, all periodic orbits having the 
same Lyapunov constant. The multiplicities in the length 
spectrum of periodic orbits grow exponentially with the 
length, like exp(Z 7 /2), where Z 7 is made dimensionless by 
referral to a scale fixed by setting the curvature to — 1; 
since action is proportional to the orbit length the action 
spectrum also has an exponentially growing degeneracy. 

The quantum energy levels for T*(2,3, 8) are found 
as eigenvalues of the Laplace-Beltrami operator. The 
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FIG. 1: Fundamental domain ofT*(2,3, 



boundary conditions can be either Dirichlet or Neumann. 
There are 2 3 = 8 quantum mechanical problems, all re- 
lated to the same classical system. In problems stem- 
ming from desymmetrization of the regular octagon, the 
boundary conditions on the triangle sides are chosen to 
obtain the spectrum for a particular irreducible repre- 
sentation. The boundary conditions on the hypotenuse 
N and the shorter leg M must then be the same, i.e., 
both Dirichlet or both Neumann; four such possibilities 
exist all of which lead to arithmetical systems with near- 
Poissonian spectral statistics. The remaining four choices 
where N and M host different boundary conditions, lead 
to pseudo-arithmetical systems with the Wigner-Dyson 
statistics of the orthogonal universality class. 

Form factor and diagonal approximation: The spectral 
form factor following from Gutzwiller's trace formula is 
a double sum over periodic orbits, 
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where S* 7 ,T 7i /i 7 , A 1 — A* are action, period, Maslov in- 
dex, and stability coefficient of the orbit 7; the Heisen- 
berg time Tjj = with A the mean level spacing, is 
used as a unit of time such that r becomes a dimension- 
less time; the angular brackets (...) denote averages over 
the energy shell and over a small r interval. 

Of special interest are the pairs of orbits obviously im- 
mune against destructive interference of their contribu- 
tions, namely with the same action and Maslov index. 
For generic chaotic dynamics these are the trivial pairs 
7' = 7 and, if time reversal is allowed, the pairs of mutu- 
ally time reversed orbits, 7' = 7™. Discarding all other 
pairs one gets Berry's diagonal approximation 0], 
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here g is the average multiplicity of the action spectrum 
which is 1 in the absence of time reversal (unitary uni- 
versality class) or 2 in the presence of time reversal (or- 
thogonal class) ; the result yields the first-order term of a 
power series in r, in agreement with RMT [J]. 



Turning to arithmetical chaos we can carry over the 
foregoing reasoning, except that the multiplicity is now 
exponentially large, g oc e'/ 2 . However, orbits with the 
same length bear no geometric similarity, and it is not at 
all obvious that their Maslov phases are the same. 

As will be shown, the Maslov phases do indeed coincide 
for all orbits in a length multiplet, a negligible fraction 
apart, in the arithmetical case. Consequently, the form 
factor exhibits almost instant increase at r > 0, similar 
to the jump of the integrable case @, K(0) — and 
K{t) = 1 for t > 0. 

In contrast, the Maslov index of pseudo-arithmetical 
systems will turn out to fluctuate randomly within each 
fixed-action multiplet. The usual diagonal approxima- 
tion with g = 2 then holds since in each multiplet only 
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terference. The diagonal approximation thus suggests 
universal spectral fluctuations. 

The Maslov index of an orbit in our billiard is de- 
termined only by the number No of reflections from 
the sides with the Dirichlet boundary condition: each 
such reflection changes the Maslov phase by tt. There- 
fore the contribution of an equal-action pair 7, 7' is 
A 2 (— 1) D D where Np, N' D respectively refer to 7 and 
7'; in long orbits both Np and N' D are large pseudo- 
random integers. We shall demonstrate that in arith- 
metical systems equal-length orbits have, in their over- 
whelming majority, No of the same parity (even or odd) 
such that contributions of all pairs of them are posi- 
tive and add up. On the contrary, pseudo-arithmetical 
systems have uncorrelated parities of Nu,N' D and the 
equal-action contributions mutually cancel, apart from 
the standard pairs of the orthogonal universality class. 

Numerical observations: The triangular billiard affords 
symbolic dynamics, and calculating its orbits involves 
generating all allowed sequences of symbols L, M, N each 
standing for the visit of the respective side. We denote 
by A 7 , [Ay, z/ 7 the number of symbols L, M, N in an orbit 
7; the total number of symbols in 7 is n 7 = A 7 + /i 7 + z/ 7 . 
Studying up to a million orbits we find: 

— Orbits within a given length multiplet A almost al- 
ways have n 7 with the same parity. We can therefore 
speak about A g - and A u -multiplets depending on par- 
ity of w 7 ,7 G A. Exceptions are extremely rare and in 
fact amount to a negligible fraction: E.g., among approx- 
imately 257000 orbits with the length £ 7 < 16 grouped 
into more than 13000 length multiplcts, only 4 multiplets 
are "git-degenerate", i.e., contain orbits with both even 
and odd number of symbols; these offenders have lengths 
I = 10.6999964, 12.2422622, 13.7571382, 15.2857092. 

— All orbits in a given A without git-degeneracy have 
A 7 of the same parity which automatically leads to def- 



inite parity of fi 7 



On the other hand, /Lt 7 and 



separately have no definite parity within A. 

These observations, in particular the rarity of multi- 
plets with git-degeneracy, suffice to explain the sensitivity 
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of the level statistics to the boundary conditions. We de- 
note by the phase jump on reflection from the side L, 
which is for the Neumann and 7r for the Dirichlet condi- 
tion on L; similarly, denote the phase jumps on 
M and N. The contribution to the diagonal form factor 
of a pair of orbits (77') belonging to the same A is 

K ' = yl2 e i(A^-Ay)*£+i(/i 7 -|*y)*M+»(f-y-I'y)*W 
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the phase proportional to disappears since, in the ab- 
sence of the git-degeneracy, A 7 — Ay is even. If the bound- 
ary conditions on M and N are different (the pseudo- 
arithmetical case) one of <&m,&n is zero and another 
one 7r. The contribution K 1 ~ j ' can then be of any sign 
and, summed over all pairs of a large multiplet, except 
the standard 7' = 7,7 TR , averages to zero. The form 
factor will be the Wigner-Dyson one for the orthogonal 
universality class, g = 2. 

On the other hand, in arithmetical systems with the 
same boundary conditions on the sides M, N we have 
<$>M = $jy = ^MNi and the contribution of every pair 
within A will be positive since p 7 + — /iy — i/y is even, 
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Eq.([2]) then applies with the abnormally high value of g 
bringing about the nearly vertical rise of K(r) at r > 0. 

Analytic reasoning: Using group-theoretical properties 
15] of T*(2, 3, 8) we have substantiated these results an- 
alytically The symbols L, M, N can be associated with 
transformations mapping the interior of the Poincare disk 
\z\ < 1 onto itself. Each elementary operation involves 
complex conjugation K : z — > z* followed by a Mobius 
transformation z — > (az + b) / (b* z + a*); the three respec- 
tive matrices ( , ) can be chosen as 
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where a, ft, 7 denote the quartic irrationalities 



V2-1, p = yV2, 
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complex plane (an invariant geodesic of the transforma- 
tion) crossing the fundamental domain; its part inside 
the domain is the piece of 7 between the two bounces 
against the sides given by the first and last symbol of 
the code. Cyclically shifting the code by one symbol one 
analogously gets the next orbit piece, and so forth flBj ] . 
The matrix p y yields the orbit length as 



2 cosh ^ 
2 

2 sinh ^ 
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Trp 7 , 
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n 7 even, 



ra 7 odd, 
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where Atr p is the sum of the off-diagonal elements of p. 

It can be shown by induction that the matrices /o 7 can 
be of two arithmetical types ([15j. Appendix B), 
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Here u k ^ Rl v k . Rl u k j,v k j are algebraic integers, 



u k , R = m k , R + n k , R V2, 
VhR = Pk,R + qk,nV2, k 



(7) 



1,2, 



with integers m k , R ,n ktR ,Pk, R ,qk,R; the imaginary parts 
have the same appearance. Appending a symbol L (pure 
complex conjugation) to the code of the orbit doesn't 
change the type of p 7 whereas appending M or N toggles 
the type, p"> «-» p^>. The matrices pu and pn belong 
to the type p^; therefore p 7 belongs to the type p^ if 
the sum /i 7 + v~ of the number of symbols M, N in the 
code of 7 is even, and p 2 if /i 7 + f-y is odd. 

Eqs. (|5|6j) entail four types of orbit lengths: Orbits 
with even numbers of symbols n 7 have lengths 



(a) 2 cosh -| = u 1:R . 



(b) 2 cosh ^ = 7 w 2 ,fl, (8) 



while orbits with odd n 7 have lengths 

(c) 2 sinh y = av hR , (d) 2 sinh y = (3 v 2 , R ■ (9) 

The integer components of u, v in these equations are 
restricted [lj| by the inequalities 



(a) \m - nV2| < 2 , (b) 
(c) \p-qs/2\<2a, (d) 



\P- 



- nV2| < V27 , 
qV2\ < V2f3 . (10) 



Let us write the code of the orbit 7 starting with an 
arbitrary symbol and multiply the associated elementary 
operators. The product is either a pure Mobius transfor- 
mation with a certain matrix p 11 if the number of sym- 
bols 7J 7 is even, or K followed by p 7 (odd n 7 ). The 
cumulative transformation leaves invariant a circle in the 



The length types (a,c) are connected with the matrices 
p-y of the type p^ whereas (b,d) are connected with p( 2 \ 
Considering the connection between the type of p 7 and 
the code of 7 we see that the type of the orbit length is 
uniquely defined by parities of the number of b ounces A^/ 
and + h> 7 , see Table [H 
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TABLE I: Orbit length types for different parities of symbol 
numbers 

A even A odd 



(/i + v) even 
(fx + v) odd 



It is elementary to prove that the numbers in r.h.s. of 
the equations (a),(b) in (j8]) cannot be equal unless they 
are zero; the same is true for (c),(d) in ((9]). Therefore a 
length multiplet can never simultaneously contain orbits 
of the types (a) and (b) or (c) and (d). On the other 
hand, an orbit with an even number of symbols n 7 can 
have the same length as one with odd n 7 , and such equal- 
ity happens for some rare combinations of Uk,R,Vk,R', 
these are the cases of git-degeneracy mentioned above. 
Equating one of l a , lb to one of Z c , Id we obtain an equa- 
tion for u, v equivalent to two diophantine equations for 
the integers m,n,p,q. (Not all solutions of these equa- 
tions correspond to really existing length multiplcts since 
(JE}, ((9]) are only necessary conditions.) A careful analy- 
sis of the latter equations (see below for details) for the 
cases l a = l c or lb = Id reveals that the solutions form 
an equidistant sequence Z 



(i) 



2arsinh2 1//4 . Similarly solutions for l a 



lb = l c are described by I 



(2) 



s k, k = 1, 2, . . . with 
Id or 
, with 



J 1 1 



1,2, 



s 11 = 2 arsinh yl + \/2. E.g., the exceptional length 
I = 10.6999964 = 7s 1 pertains to two multiplets, (b) 
with m = 138, n = 97, and (d) with p = 88, q = 63. 

Since the exceptional lengths appear in equidistant se- 
quences the multiplets with git-degeneracy are exponen- 
tially outnumbered, as the length I grows, by the multi- 
plets of definite type and hence parity of A and fi + v. 
The git-degenerate multiplets can thus be discarded for 
lengths corresponding to a finite fraction of Tjj . 

To conclude, we have shown that "all" (all save for 
a negligible fraction of multiplets) orbits with the same 
length of the T*(2, 3, 8) billiard have the same parity of 
the number of bounces against the longer leg L of the 
triangle; this is also true for the parity of the sum of the 
number of bounces against the sides M, N. As a result 
"all" orbits with the same length in the arithmetical case 
have the same Maslov index such that all diagonal terms 
in the form factor are positive. A similar mechanism 
must exist in all other arithmetical systems with non- 
vanishing Maslov phases. In pseudo-arithmetical systems 
the Maslov phases of degenerate orbits are uncorrelated. 

Off- diagonal corrections: We based our reasoning on 
the diagonal approximation. But the essence of our ar- 
gument carries over to the higher-order terms of the r- 
expansion valid for r < 1 0] as well as for the behavior 
at times exceeding the Heisenberg time as far as the 
pseudo-arithmetical case is concerned: The high multi- 
plicity of length multiplets of orbits is rendered irrelevant 
by destructive interference of random Maslov phases. 



For arithmetical systems the interplay of high length 
degeneracy and bunches of orbits that very nearly coin- 
cide in configuration space, apart from reconnections in 
close self-encounters, may be more difficult to capture. 
On the other hand, the Heckc symmetries intuitively sug- 
gest (near) Poissonian level statistics. 

Technical note on exceptional multiplets: We briefly 
indicate how existence and uniqueness of the exceptional 
lengths can be ascertained, for the four possible cases 
la = Ic-Ja = Id, h = l c ,h = Id which we shall refer to 
as ac, ad, be, bd. Starting with the case ac we infer from 
Eqs. (|8I9|) that the exceptional lengths obey 



av 



av 



(11) 



Now take two solutions uq + avo and u + av correspond- 



ing to the lengths l ,l. Then since exp [± 
exp (±y) ex P (^l)' the product 



uq + avc, it + av 



>)] 



(12) 



is a solution with length I' = Z + Z , provided the integers 
in the rhs are divisible by 2. 

The numerically found solution of smallest length Zq = 
2s 1 = 3.05714183896200 was u a = 2 + 2^2, v = 4+2^2, 
or mo = 2, no = 2,po = 4, qo = 2. Substituting the latter 
into (11211 we obtain 



m! = m+ 2n + 2q, (13) 

n = m + n + p, 

p' = 2m + 2n+p + 2q, 

q' = m + 2n + p + q; 

obviously if m,n,p,q are even then so are the primed 
numbers. We thus face the sequence of solutions 
( " 0+ 2 at, ° ) fc with the equidistant lengths kf Q . 

The transformation (fT3"|) can be inverted; the doubly 
primed numbers 



n — 
ii 

P = 

q" = 



m + 2n — 2q, 
m + n — p, 
—2m — 2n + p + 2q, 
—m — 2n + p + q: 



(14) 



yield an orbit of length Z" = Z — Zq such that a ladder of 
decreasing lengths is obtained. 

In order to make sure that the transformations (jl3U4[) 
really yield orbits we must check that the inequalities 
(a,c) in (ITOl) are not violated. To that end we note that 
the variables x = (m — nV2)/2, y = (p — qy2)/2a span 
an invariant subspace of the transformations (|13|14p ; the 
ensuing transformation (x,y) — > (x',y') is a two dimen- 
sional rotation preserving x 2 + y 2 . Since the aforemen- 
tioned minimal-length solution mo = 2, no = 2,po — 
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4, q = 2 has the property x 2 +y 2 = 1 and since that "nor- 
malization" is preserved we indeed conclude the preser- 
vation of the inequalities under study. 

Finally, let us demonstrate that there are no lengths 
of the type ac outside the equidistant-length ladder just 
established. Momentarily assuming the existence of such 
a freak length we can apply the length reducing transfor- 
mation (|14[) until arriving at a reduced length I" £ [0, Iq] 
which must obey 

, I" m" + n"V2 i= u Zo 

1 < cosh— = — < 1 + V2 = cosh— , (15) 

2 2 2 K ' 

< sinh - = (p" + vW) - < (2 + V2) a = sinh | . 

On the other hand, the conservation of x 2 + y 2 sub- 
jects the freak to (x"f + [y"f = (m" - n" ^ /A + 
(p" - q"V2) 2 /4a 2 = x 2 + y 2 < 2, and thus to 

1 to" - n"V2\ < 2V2, \p" - q"V2\ < 2^2a. Upon 
checking the remaining finite number of quadruplets of 
integers m" ,n" ,p" ,q" with lengths inside the interval 
< I" < Iq we conclude that no freak can exist. 

The discussion of case ad proceeds in close analogy 
with the previously treated ac. The only differences are 
the replacements of (i) the quartic irrationality a by j3 
and (ii) the minimal length by l 1 ^ = 2s 11 = 4.8969, 
the latter corresponding to uq — 6+4V2, «o = 4+4^2- 
All other lengths are given by Hq 1 , k = 1,2, ... . The 
proof repeats the previous case word for word. 

To treat case bd we employ 7 = aj3 to write the equa- 
tions for the lengths as 

1 (3 (au + v) _i /3(au — v) , < „. 
e2 =^^ ~> e2= ^ ■ ( 16 ^ 

The smallest length corresponds to u m — 2 + \[2, v m = 
\]~2 and is exactly s 1 , i.e., half of the smallest length of 
case ac. Solutions of (fTl)|) do not have in general the group 
property of the case ac; in particular, if I is an exceptional 
length of the type bd then 21 is not. Indeed, squaring the 
rhs of (HU we get y/2( au ± v) /4 which cannot be an 
algebraic number of the type j3 (au + v) /2. On the other 
hand, all odd multiples of I do belong to the admissible 



type. The same reasoning as in the cases ac, ad shows 
that all solutions can be represented as (k + 1/2) Iq, k — 
0,1,.... 

The remaining case be is related to bd by the replace- 
ments a ■(-> (3 and Iq — > Iq 1 . All solutions can be written 
as (fc + 1/2)^, k = 0, 1,.... 
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